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Abstract: 

ON 

in 

We construct an approximate expression for the cross section for Higgs production in gluon fusion 
at next-to-next-to-next-to-leading order (N 3 LO) in a s with finite top mass. We argue that an accurate 
approximation can be constructed by exploiting the analiticity of the Mellin space cross section, and 
the information on its singularity structure coming from large N (soft gluon, Sudakov) and small N 
(high energy, BFKL) all order resummation. We support our argument with an explicit comparison 
of the approximate and the exact expressions up to the highest (NNLO) order at which the latter are 
available. We find that the approximate N 3 LO result amounts to a correction of 17% to the NNLO 
QCD cross section for production of a 125 GeV Higgs at the LHC (8 TeV), larger than previously 
estimated, and it significantly reduces the scale dependence of the NNLO result. 
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1 Introduction 



The dominant Higgs production mechanism at the LHC is gluon fusion via a heavy fermion loop 
(mainly a top quark) and indeed the recent announcement of the discovery of a Higgs- like parti- 
cle [2II3] is largely based on events in this channel. In view of this, an accurate determination of the 
cross section in this channel is of great interest. Next-to-leading order (NLO) corrections to the inclu- 
sive cross section, originally computed in Refs. [4j[5] in the large top mass (mt — > 00) approximation, 
and in Ref. [6] for general m< are known to be as large as the leading order, and the NNLO corrections 
(first computed in Refs. [7]-[9] in the mt — > 00 limit and for finite top mass in Refs. |10f[15] ) about 
half as large as the leading order. The significant scale dependence of the NNLO result suggests that 
corrections at yet higher orders are not negligible: in fact they currently account for half or more of 
the uncertainty on the theory prediction for the cross section [IB] (the other half being due to parton 
distributions and the strong coupling). 

While computations of the full N 3 LO correction to the cross section are in progress 
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it 



is interesting to derive approximate expressions for it. Several of us have argued (see e.g. 20-22]) 
that accurate approximations to partonic cross sections may be obtained from knowledge of their N 
space singularity structure, both at finite perturbative order, and at the resummed level. Because the 
N — > 00 singularity and the rightmost singularity at finite N are known to all orders in a s respectively 
from threshold (Sudakov) and high energy (BFKL) resummation, if this is indeed the case it is possible 
to construct reliable approximations even to very high orders in a s . The possibility of constructing 
approximations based on the combination of results from large and small N resummation has also 
been considered in 123,241. 



In this paper, we will pursue this idea in the context of Higgs production in gluon fusion: we will 
determine the dominant small N and large N singularities up to N 3 LO from resummation arguments, 
and, after testing our methodology against known results up to NNLO, we will use them to construct 
a N 3 LO approximation. 



2 The partonic cross section and its singularities 

The factorized Higgs production cross section is 

. V-^ f 1 dZ (T n\ 1 A / o . n. m 2 rr m 2 rr\ uAr . ^ 

^ J T Z \Z J Z \ flj^J s 

where Jfij(z,fi 2 ) are the parton luminosities 

^(z,n 2 )= f 1 —fi(-, M 2 )/,-(x,/i 2 ). (2-2) 
We introduce coefficient functions Cy, defined as 

&ij(z,m 2 H ,a a (nD, T ^f, V ^f \ = zao(m 2 H ,a s (vl)) Cij(z,a s (p%), (2.3) 
V Mf Mr / V Mf Mr / 

where o"o is the leading order (LO) partonic cross section, so that the coefficient function is normalized 
to 6(1 — z) at leading order: 

Cij(z, a.) = 5(1 - z)5 ig 5 jg + a s c\f(z) + a 2 s cg\z) + a 3 cf(z) + <D(a% (2.4) 

and for simplicity, we have suppressed the dependence on renormalization and factorization scales 
/i F ,/U R . In the sequel, we will concentrate on the gluon fusion subprocess, while the contribution from 
other subprocesses will be only briefly discussed in Section [4j so in most of the discussion below we 
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will drop the parton indices ij, and assume that both the coefficient function and luminosity refer to 
the gluon channel. 



Because the cross section Eq. (2.1) is a convolution, its Mellin transform 

a(N,m 2 H )= [ dTT N - 2 a(T,m 2 H ) (2.5) 
Jo 

factorizes in terms of the Mellin space luminosity and coefficient function, respectively defined as 

JSf(JV) = f dzz N - 1 ^(z) 
Jo 



C(N,a s ) = I dzz N - 1 C{z,a s ), (2.6) 

according to 



o 



a(N,m 2 H ) = a {m 2 H ,a s )J?(N)C(N,a s ). (2.7) 

While in momentum space the coefficient functions are distributions, if the Mellin transform in- 
tegral has a finite convergence abscissa, the N space coefficient function is an analytic function of 



the complex variable N, given by the integral representation Eq. (2.6) to the right of the convergence 
abscissa, and by analytic continuation elsewhere. Therefore, it is fully determined by knowledge of its 
singularities. 

The singularity structure of the perturbative expansion of C(N,a s ) is relatively simple. At any 
perturbative order, the rightmost singularity is a multiple pole located at N = 1 [25], with further 
multiple poles along the real axis at N = 0,-1,-2,..., with residues of order one (this is also what 
is found in all known fixed order calculations); Re TV = 1 is the convergence abscissa of the Mellin 
transform, and as N — > oo, C(N, a s ) grows as a power of In A. While knowledge of the residues of all 
poles is required in order to fully determine the function C(N,a s ), its behavior in the physical region 
1 < Re A < oo is mostly controlled by the residues of the leading (rightmost) pole at A = 1, together 
with that of the singularity at infinity. Both are known from resummation: Sudakov (soft gluon) 
resummation determines to all orders in the strong coupling the coefficients of the ln m N terms which 
control the behavior as iV — > oo, while BFKL (high energy) resummation determines the residues of 
the leading ^ N ^„ multiple poles. 



This suggests that an approximation of the coefficient function Eq. (2.6) may be constructed by 
simply combining the large iV (soft) and small A (high energy) terms, 

C appr0 x(A, a„) = C soft (A, a„) + C h . e .(N, a s ), (2.8) 

where C so f t contains terms predicted by Sudakov resummation and Ch. e . terms predicted by BFKL 
resummation. It is clear, however, that this is only correct if the small A singularities, controlled 
by Ch.e., are unaffected by C so ft, while the large N logarithms, controlled by C so f t , are unaffected 
by Ch.e.- This is clearly nontrivial: for example, a term proportional to ln m A has a cut at A = 0, 
while at each fixed order the expected behavior of the coefficient function is a pole, rather that a cut. 
So the approximate expressions for C so f t and Ch. e . should reproduce this behavior, with no spurious 
singularities. 



We will show in the sequel that an approximate expression of the form of Eq. (2.8) is possible, but 
both C so ft and Ch. c . will have to be carefully constructed. Indeed we will now show that constructing 
C so ft in such a way that the small A singularity structure is preserved, the agreement at large A is 
considerably improved. This result may seem surprising, but it is in fact a consequence of analyticity. 

2.1 Large N 

We first discuss the computation of the large A (soft) part of the coefficient function. All contributions 
to C(N,a s ) which do not vanish as A — > oo may be computed from Sudakov resummation, using 



3 



techniques summarized long ago in Ref. [26]. The resummed coefficient function has the form 

1 



C res (N, a s ) = go{a s ) exp 



- gi (a s In N) + g 2 (a s In N) + a s g 3 (a s hxN) + ... 



with 



go(a s ) = 1 + a s g ,i + a 2 s g 0t2 + 0(a 3 s ), 

oo 
k=k oi 



H),i>2 



(2.9) 

(2.10) 
(2.11) 



Inclusion of all gi with 1 < i < + 1 and of go U P to order a k_1 gives the next fc -to- leading log 
approximation to In C res (N,a s ); it determines the coefficient of all contributions to the coefficient 
function of the form a™ ln m N with 2(n— k) + l < m < 2n. This can be extended to 2{n — k) < m < 2n 
by also including the order a k contribution to go . The functions g\ , g 2 and g% are known exactly, while 
go is known up to 0(a 2 ). The function g± is only known in part 27,28 , but the missing information 
(the 4-loop cusp anomalous dimension) only enters at O(af). We can thus determine all large N 
non- vanishing contributions to C(N,a s ) up to 0(a 2 ), and all logarithmically enhanced contributions 
(but not the constant) to 0{a 3 ). 

The accuracy of an approximation to the Higgs production cross section at the LHC based on 
the dominance of threshold terms can be studied |29| by using the saddle point method to determine 
which is the region in N space that gives the bulk of the contribution to the cross section. It turns out 
that, despite the fact that Higgs production at the LHC is far from the kinematic threshold, partly 
because of the underlying partonic kinematics and partly because of the shape of the cross section, at 
the LHC with 8 TeV center-of-mass energy, logarithmically enhanced terms are still providing most 
of the cross section, though the situation gradually changes as the center-of-mass energy increases. 

However, our goal here is to construct an approximation to the coefficient function which holds 
for all N in the physical region. Now, it has been observed long ago [30| that the the quality of the 
soft approximation to the full coefficient function significantly depends on the choice of subleading 
terms which are included in the resummed result: indeed, while resummation uniquely determines 
the coefficients of logarithmically enhanced terms, there is a certain latitude in defining how the soft 
approximation is constructed, by making choices which differ by terms which vanish as N — > 00. 
A similar situation has been observed recently in Drell-Yan production at the LHC (2l], for which 
the threshold approximation is generally expected to be less good than for Higgs production. By 
comparing results which differ by terms which vanish as iV — > 00, we will now show that several 
preferred choices for such subleading terms are favored by the requirement that some aspects of the 
known small iV singularity structure of the exact result be reproduced. 

In order to outline our strategy, let us work with the simplest example. Let us first suppose that 
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we know the N space resummed coefficient function and that we want to extract from Eq. (2.9) an 
approximate expression for the 0(a s ) coefficient C^\z), which is given by [6,; 



C«(z) = 4A g (z)V 1 (z) + dd(l -z)- 2A g {z)^- + K gg (z), 

V(i-*)\ 



V k (z) = 



1-z / ' 
\ / + 

C A 1 - 2z + 3z 2 - 2z 3 + z 4 



7T 



(2.12) 
(2.13) 

(2.14) 



The constant d and the function lZ gg (z) are known functions of m^/mt; in particular 1Z gg (z) is an 
ordinary function, regular in z = 1, so its Mellin transform vanishes as iV — > 00 and therefore its 
specific form is of no relevance for the large N behavior. 
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Expanding Eq. (2.9) to 0(a s ), and keeping NLL terms, we find 

C ies (N, a s ) = l + a s C£}(N) + 0(a 2 s ), (2.15) 
C«(Af) = gi, 2 In 2 N + 5 2,i In AT + <?o,l, (2.16) 

with 

2C A AC A 
5i,2 = , 92,1 = 7s, ( 2 -!7) 

where 7 B is the Euler-Mascheroni constant. The asymptotic behavior of the 0(a s ) coefficient as 
N — > oo is correctly reproduced by this expression, in that 



lim 

TV-s-oo 



0; 



(2.18) 



where C^(N) is the Mellin transform of Eq. (2.12); the constant <?o,i is fixed by this condition 



On the other hand, the behavior of Eq. (2.16) at small values of N is incompatible with the known 
singularity structure. In particular, there is a logarithmic branch cut starting at N = which is 
definitely unphysical, as the exact coefficient function has poles and not cuts at small N. This cut is a 
subleading singularity, given that the leading singularity is located at N = 1, but close enough to the 
leading one that the behavior of the coefficient function can be significantly affected. Even if we plan 
to eventually improve this expression by introducing the correct singularity at JV = 1 according to 



Eq. (2.8), the logarithmic singularity will interfere with it and spoil the accuracy of the approximation. 



This problem, however, is an artifact of the large N approximation, since powers of In N are the 
large N approximation of powers of the digamma function tpo(N) appearing in fixed order computa- 



tions. Indeed, the inverse Mellin transform of Eq. (2.16) (using Eq. (A. 6b) of Appendix A.l) is seen 
to be 



C«(z, a s ) = g Q ,iS(l - z) + 2g 1)2 V l ° s (z) + (2 7b51i2 - g 2)1 )V l ° s (z), 



7T 



where 



v l ?Hz) 



1 



ln K In 



In 



i 



(2.19) 
(2.20) 



which is seen to differ from the soft contribution Eq. (2.13) to the exact result Eq. (2.12) 



This can be understood noting that singular terms as z — > 1 arise from integration of the real 
emission diagrams over the transverse momentum of the gluon, which has the form 



(2.21) 



where A is a collinear cut-off and p gg (z) is the LO gluon-gluon Altarelli-Parisi splitting function for 

z < 1, 

A 9 {Z) 



Pgg( z ) 



1 - z 



(2.22) 



with A g (z) given by Eq. (2.14). 



Indeed, Eq. (2.21) shows that logarithmically enhanced soft terms, rather than being proportional 
i 



to — 7— t- 5 - , are of the form 

m ± 



1 Ini^ 



1 - z 



1-z I 



ln(l - z) + 0(1 - 2) 



(2.23) 



5 



and they appear with a coefficient proportional to the Altarelli-Parisi splitting function. Explicitly, 
the latter in the z — >• 1 limit may be expanded as 



A 9 (Z) 



l-(l-z) + 2(l-z) 2 + 0[(l 



(2.24) 



C A 

^z) = — 

Logarithmically enhanced contributions to the coefficient function are generated by the first terms in 
both expansions Eqs. (2.23) and (2.24), namely ln(l — z) and A g (l) respectively. 



We will now argue that an optimal choice of the soft approximation, differing from Eq. (2.19) by 
subleading terms, is obtained by writing the large soft logs as powers of In , so in particular retaining 

the y/z in the denominator despite the fact that it is subleading, and furthermore, by retaining at 
least the first correction on the right-hand side of Eq. (2.24), also subleading. Therefore in this case 
our suggestion consists in the simple replacement 

A g (l)V^(z)^A g ^ m (z)V l (z) 



(2.25) 



in Eq. (2.19), where A g ^ m (z) is a finite m-th order expansion of A g (z) about z = 1, Eq. (2.24), and 

ln(l — z)\ In 



V x {z) = 



1 



1 



(2.26) 



Note that we have chosen to apply the plus prescription only to the first term, singular in z = 1, which 
is the natural choice in fixed order calculations. In this way, T>i(N) differs from T>i(N) only by terms 
vanishing at large N. Since T>° S [{N) and T>i(N) differ at large N by a constant, the coefficient 50,1 
must be modified accordingly, in order that the requirement Eq. (2.18) be satisfied. These technical 



details are discussed in Appendix A.l 



Our conclusion Eq. (2.25) relies on the following arguments: 
• The replacement of T>^° s (z), whose Mellin transform is 

v\° s (N) = l[ln 2 N + 2 lE lnN] 



(2.27) 



with t>x(z), whose Mellin transform is 

1 

2 



A(N) = - [^(N) + 2 7b Vo(A0 + C2 + r 



(2.28) 



removes the logarithmic branch cut of T>° g (N), which is incompatible with the known analytic 
structure of the coefficient function. The only singularities are now isolated poles, as in the exact 
expression. 



• The same features are shared by the Mellin transform of T>\(z), that is 

Vi(N) = \ [i,l{N) - MN) + 2 lE MN) + C2 + T 



(2.29) 



However, the presence of ipx(N) exactly cancels the double poles of iPq(N) in N = 0, — 1, —2, . . ., 
which are there in the exact result. Therefore, the choice of T>\{N) is preferred over T>i(N). 



In the replacement Eq. (2.25) the factor A g (z) is expanded up to a finite order m > about 
2 = 1. This is because the inclusion of the full A g (z) would introduce a spurious singularity in 
N = 1. Indeed, the Mellin transform of A g {z) ~Di(z) is given by 



dz z N - 1 A g {z)V 1 {z) 



Ca 

7T 



T>i(N - 1) - ZDi(N) + 3£>i(Ar + 1) - 2Pi(A^ + 2) + T> X (N + 3) . (2.30) 
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The first term, due to 1/z in A g (z), has a double and a simple pole in N = 1, while the 
exact singularity is a simple pole, with a (ma /nit )-dependent coefficient controlled by small z 
resummation. The expansion of A g (z) in powers of 1 — z to any finite order is not singular in 
z = 0, and therefore does not affect the singularity structure around N = 1. 



We turn now to the general case. Each of the above arguments can be generalized to all orders, 
where N space resummed results contain powe rs of ln fc N, whose inverse Mellin transform is a linear 



combination of distributions U° g (z) Eq. (2.20) with j < k — 1. The fact that the NLO result in z 
space depends on powers of In rather than In In - is of kinematical origin and ultimately comes 



from the upper bound for the transverse momentum of emitted gluons, Eq. (2.21), and therefore it 



persists to all orders. It follows that the exact result to all orders is expressed in terms of distributions 
T>k{z), defined in Eq. (A. 2c) of Appendix A.l in analogy with Eq. (2.26). The Mellin transform of such 



distributions, T> k (N), first, has poles rather than cuts as small N singularities, and also, in comparison 
to the distributions T>k(N), lacks contributions proportional to powers of tp k (N) with k odd, which 



would change the pole structure (see Appendix A.l). 
It has been shown in Refs 
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that the factor A g (z) 



Eq. (2.14), is present to all orders, 



because the full leading order anomalous dimension exponentiates. However, terms beyond the first 
in its expansion Eq. (2.24) generate contributions a™(l — z)i ln 2n ~ 1 (l — z) with j > to the coefficient 
functions, which are generally of the same order as other terms which we do not control. However, 
it can be shown [33] that the inclusion of the 0\{1 — z) 1 ] term in the expansion Eq. ( 2.24| ) correctly 
predicts, after exponentiation, the subdominant contributions of the form a"ln 2n_1 (l — z) (i.e., in N 
space, terms behaving as a™N~ l ln 2n_1 N at large N) to all orders, so the inclusion of this term rests 
on firm ground. 



Including the 0[(1 - z) 1 ] from Eq. ( |2.24[ ) we get 

C A 



A 9,l( Z ) 



7T 



[l-(l-z)] = zA g (l), 



which is easily implemented to all orders by the replacement 

v[°Hz) -4 zV k (z); V[°HN) -4 V k (N + 1). 

Including also the next order gives 



which amounts to replacing 



C A 



IT 



[1 - (1 - z) + 2(1 - z) 2 } = [2 - 3z + 2z 2 ]A g (l), 



2V k {N) - 3V k (N + 1) + 2V k (N + 2) 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



in the iV space expressions. The third order term of the expansion of A g (z) is accidentally zero, so 
A g ^{z) = A g ^(z). We have checked that the inclusion of terms of order (1 — z) 4 and higher in the 
expansion of A g {z) does not affect our results significantly. We will consider both the expansions 
to first and second order, and use their difference as a means to estimate the uncertainty on the 
result. Specifically, we will take the mid-point between them as our best prediction, with the first- 
and second-order expansion result giving the edges of the uncertainty band. 

In summary, our soft approximation (to be combined with small iV terms determined in the next 



Section) is constructed in the following way. The resummed expression Eq. (2.9) can be rewritten 



C rcs (N,a s ) = g (a s ) exp ]T a?£ b n>k ^(iV), 



(2.35) 



n=l 



k=0 



where the coefficients b U:k are obtained from the functions gi, Eq. (2.11), and have been determined 
up to n = 3 28 . The function go(a s ) is known only up to 0(a 2 ); the uncertainty associated to 50,3 
will be discussed in Sect. 
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The replacements Eq. (2.32) or (2.34) are then applied to Eq. (2.35). We obtain, respectively, 

(2.36a) 



Csofti (N, a s ) = go(a s ) exp 



n=l k=0 

oo n 



Csoft 2 (N, a s ) = g (a s ) exp ^ a™ ^ 

n=l k=0 



2V k (N) - 3V k (N + 1) + 2V k (N + 2) 



(2.36b) 



where we have defined 



g (a s ) = g (a s ) exp 
d k = lim 



7V-s>oo 



^2 a ™ bn > kdk 

n=l k=0 

V k (N)-V k °Z(N)} 



(2.37) 
(2.38) 



so that the condition Eq. (2.18) is satisfied to all orders after the replacement. Explicit expressions 



for the coefficients b n>k and d k are given in Appendix |A.1[ 



Equations (2.36) can be cast in the form 

oo 

C soit (N,a s ) =g {a s ) exp ^ a™ S n (N), 



n=l 



which is now expanded in powers of a s : 

C soit (N, a s ) = l + a s cV( N ) + a 2 s C^ t (N) + a? t C®(N) + 0(a A s ) 



We obtain 

c£l(N) = S 1 (N)+g 0A 

C S oftW = ^?(^) + S 2 (N) + go,iSi(N) + g ,2 



(2.39) 

(2.40) 

(2.41a) 
(2.41b) 



C SU N ) = l S f( N ) + Si(N)S 2 (N) + S 3 (N) + 9 o,i (^SfiN) + S 2 (N)\ + go^N) + g , 3 . (2.41c) 



As a test of our procedure we now compare the first two orders of our soft approximations Eq. ( 2.36 ) 



to the full result. Note that in the sequel when comparing to known results, and also when constructing 
our O(a^) approximation, we will always be retaining the exact mt dependence. 

As terms of comparison, at NLO we use the finite-m^ result of Ref. j6j (using the numerical 
implementation of Ref. |31]), while at NNLO we use the approximate finite-mt result obtained by 



matching the double expansion in powers of 1 — z and mH/mt of Refs. 11 12 to the known small z 
terms computed in Ref. [To] according to Ref. [13] (see Refs. [14, 15 for further approximate finite- mt 
results). Note that the soft limit only depends on mt through the function go(a s ) of Eq. (2.9). 



Results are shown, as functions of N along the real N axis, in Fig. [1} We find the comparison 
in N space to be most instructive, because the coefficient function is then an ordinary function, 
rather than a distribution as in z space. Furthermore, the saddle point which dominates the Mellin 
inversion is on the real axis [29] . All this said, it should be kept in mind that the physical cross section 
is obtained by Mellin inversion of the product of the space coefficient function and luminosity: 
therefore, agreement on the real axis is certainly necessary, but in general not sufficient for agreement 
of the physical results. In particular, spurious singularities (and in particular spurious cuts) may 
substantially modify the behavior of the coefficient function in the complex plane. 

In order to understand the role of various subleading terms, we also show in Fig. [l] the results 
obtained expanding the resummed expression, Eq. (2.35), which is built up from the distributions 



S 




Figure 1. The partonic coefficient function Eq. (2.40) in N space for ran = 125 GeV at NLO (left) and 
NNLO (right) and its soft approximation in various forms: our preferred choices Eqs. (2.36), denoted as 
softi and soffo, and the simpler approximations based on T> lo6 (N) as in Eq. (2.35), denoted as JV-soft, and 
based on V(N), as in Eq. (2.42), denoted as soft-0. 



2?! og (z), Eq. (2.20), and thus it has spurious cuts starting at N = (labelled as iV-soft), and the one 



obtained expanding the resummed expression Eq. (2.35) in powers of a s and then replacing 



T%*{N) ->• V k (N) 



(2.42) 



(and adjusting the constant term), where T>k{N) are the Mellin transforms of the T)j t (z) distributions, 



Eq. (2.13), so it does not include the contributions coming from the in the phase space integration 



(labelled as soft-0). 

While our preferred options clearly provide the best approximation to the exact result in the 
soft region N > 2, it is interesting to observe that the iV-soft form, based on Eq. (2.35), despite 



having the wrong singularity structure (cuts rather than poles at small N), still provides a reasonable 
approximation, though not quite as good as our preferred ones. This can be understood noting 



that 21 



ln fc In 

^nT 



In 



fcl 



i + o[(i 



(2.43) 



which means that T*^ g (z) 
logarithms is very similar to our C so fti approximation Eq. (2.36a), up to a shift in N by 1/2 



y/zVk{z) + delta terms + 0[(1 — z) ln fc (l — z)\ , i.e., this choice for the 



2.2 Small N 



The leading small N singularities for the Higgs inclusive cross section have been determined to all 
orders in a s in Ref. |34| in the mt — > oo limit, and in Ref. [To] for finite mj. These results have been 
obtained by means of the so-called high energy or kt factorization technique of Ref. [35] , which has been 
subsequently used to compute high energy cross section for an increasing number of processes p6[|40] 
and more recently extended to rapidity distributions (and also used to determine all order results for 



Higgs production) in Ref. 41 



In this formalism, small N singularities are obtained to all orders by computing the leading order 
partonic cross section for the relevant process, but with off-shell incoming gluons. They are extracted 
from the off-shell coefficient function C ff- S heii, defined through 

Vg*g*^H = 2 00 C ff-siiell(2, (1,^2)1 (2.44) 

where z is the scaling variable defined in the previous section, while £j = \k ti \ 2 /m'jj in terms of the 
off-shellness kf = — |&ti| 2 of the i-th incoming gluon, and the angle between the incoming transverse 
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momenta is integrated over. To do this, one defines the impact factor h(N, Mi, M2) according to 



h(N, Mi, M 2 ) = M x M 2 R(Mi) R(M 2 



■m 



2 x Mi+M 2 
H 
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I "T - ?1 / ~^?2 ° off-shell {Z, £l, £2 J 

Jo a Jo « 



dz 

x / — z~ 
lo * 

00 

J2 <HiM (mum H , Mf) M^M* 2 + 0{N - 1), 
ii,i2=0 



(2.45) 



where the pre-factor R accounts for factorization scheme dependence 42 , and in MS is given by 



R(M) = 1 + -C3M 3 + 0(M 4 



(2.46) 



The determination of the coefficients Ci lt i 2 (m,t, mn, Hf) has been reduced to quadratures to all orders 



in Ref. 10 ; they have been numerically determined up to and including second order in a s in [10] and 



up to and including fourth order in 43 



The leading singularities of the partonic coefficient function are obtained by identifying the Mellin 
variables Mj with the anomalous dimension 7+. This, in turn, is the eigenvalue of the singlet anomalous 
dimension matrix which contains, to all orders in a s , the contributions with the highest powers of the 
rightmost N space singularities. Indeed, as well known, only one of the two eigenvalues (which 
henceforth we will refer to as the "large" eigenvalue) has singularities at N = lQ while the other has 
singularities at N = 0. 



In other words, the leading singularities are found letting Mi 



Ts 



with 



= 1 



a s (m 2 



Hi 



N — I 



(2.47) 



where the coefficients e n -i^- n are determined 44 using duality [45] from the leading order BFKL 
kernel. The first 35 coefficients e r 



-■n, —n 



are tabulated in Ref. 46 ; the first few have accidental zeros, 
and are given by en,-i = C^/V, 61,-2 = ^2,-3 = 0, 63,-4 = 2C3(C' y 4/7r) 4 , 64,-5 = 0. It follows that 
to fc-th order in a s the coefficient function has a fc-th order pole in N = 1. Note that in the heavy 
top limit the small iV singularity structure is different, in that at each extra order in a s the order of 
the pole increases by two units |34| . However, these double poles are unphysical, and follow from a 
breakdown of the large mt approximation at high energy: we will thus not discuss them further. 



It has been shown in Refs. 20 47 that the nature of the small N singularity of coefficient functions 



at the resummed level is entirely determined by the singularity of the resummed anomalous dimension 
7 + . However, reproducing the correct all order small N singularity of the anomalous dimension (which 
is a simple pole to the right of TV = 1, but close to it) requires [48] the all order inclusion of two classes 



of subleading terms on top of the leading (or next-to-leading) singularities Eq. (2.47): namely, running 
coupling corrections, without which the small N leading singularity would be a square-root cut instead 
of a simple pole |49|, and anticollinear terms [50] without which the perturbative expansion of both 
the position and residue of the above simple pole would not be stable (similar conclusions can also be 



arrived at from a study [51 52 of the BFKL 25 equation) . The inclusion of a further series of all 
order running coupling corrections in the coefficient function is further required [20[|47| in order for this 
not to develop extra spurious singularities. When expanded out in perturbation theory, these running 
coupling corrections correspond to series of contributions of increasingly low logarithmic order (i.e. 
increasingly subleading): we will retain both up to the NLL order, i.e. keeping not only the leading 
singular contribution to each order in a s , but also the first subleading correction, i.e. to order both 
the contributions with a k-th and a (k — l)-th order pole in N = 1. 



1 Note that in the small iV literature, and specifically in Refs. |l0[ 43 the variable iV is usually shifted by one unit, 
so that the singularities of j + are located at N = 0, by taking x instead of x"- 1 as a kernel of the Mellin transform 



Eq. (2.61. Throughout this paper we adopt instead the more common convention of Eq. (2.61 
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For anomalous dimensions this is simply done by including the full next-to-leading singular con- 
tributi on to them. For coefficient functions, these running co upling corrections are found by letting, 



in Eq. (2.45), M. 



, with 



k+l 



7r, 



given recursively by 



20 



47 



7res I 1 + 2 
\ 7res , 



(2.48) 



where 



Ives = -/3oas^7rt SJ ( 2 - 49 ) 

and with 7+ s we have denoted a form of the large eigenvalue which includes at least the leading 



singularities Eq. (2.47), but may include other subleading contributions. 

We can now compute the small N approximation to the coefficient function. We expand the 
anomalous dimension to fixed perturbative order 



7 + = as7 (0) + ah (1) + ah (2) + O(a!). 
The leading and next-to-leading singularities of the anomalous dimension 7+ are given by 

7 ( 



7 



(i) 



N —1 

ei-2 



(N - l) 2 N — 1 
e 2 ,-3_ e 2 ,-2 Q((N _ ir n 
1)3 + (AT _l)2 + )' 



(N 



(2.50) 

(2.51a) 
(2.51b) 
(2.51c) 



where the coefficient of the leading poles can be read off Eq. (2.47): eo,-i = and e\-2 = e 2j _3 
The other coefficients are: 



0. 



eo,o 



ei - 



-llC A + 2n f (2C F /C A -l) 
12vr 

/ 13Cp 23C A \ 



18vr 2 



36vr 2 



n 



e2,-2 



C A ( 3 , UCl( 2 395C 3 ( (C\C2 71C 2 A C F C A ( 2 , 71C F C A \ 



+ 



+ 



+ 



n f . 



(2.52a) 
(2.52b) 

(2.52c) 



2vr 3 12vr 3 108vr 3 V 6tt 3 108tt 3 3tt 3 54vr 3 
The N — > 1 result for the partonic coefficient function in the gluon channel can be then obtained 



by substituting Eq. (2.50) into Eq. (2.45): 



C A bf(N, a & 



71=1 







+ «2 


7 




7 



ii,i 2 >0 

q s 2ci,o7 (0) 



+ a\ 



(2c 2 ,o + ci,i) 7 (°) 2 - 2c 2 , /3o7 (0) + 2ci, 7 (1) 
(C3,0 + C2,l )2 7 (0) ' - (3c 3 ,0 + C2,l )2/3 7 (0) +4c 3 ,o/3 2 7 



(0) 



+ (2c 2 ,o + cii^Pty 1 ) - 4c 2i0 /3o7 (1) + 2c lj0 7 (2) 



+ 0(at), 



(2.53) 



where we have omitted the dependence of the coefficients on nit, fnii and /u F for simplicity. 
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Figure 2. Comparison of the NLO (left) and NNLO (right) exact coefficient function both for finite m t 
(exact) and in the m t — > oo limit (large mt), and several small N approximations to it: the NLO and NNLO 
contributions to C A bf Eq. flgjg} (ABF), to C AB F-sub Eq. ( pEsIfr (ABF-sub) and to C h . e , Eq. p755] ) (high 
energy). 



Because we wish to combine the small N behavior which we are determining here with the large A" 



behavior determined in Sect. 2.1 we must make sure that the small A" contribution to the coefficient 
function vanishes as N — > oo. However, the coefficient function Cabf(A') Eq. ( 2.53 ) manifestly does 
not vanish in the large N limit, because of the constant contribution to 7^°- ) Eq. ( 2.51 ) which propagates 



into Cabf(A^) to all orders in a s . 

Therefore, we construct an improved small A" approximation to the coefficient function as a sub- 
tracted version of Cabf(A'). The subtracted coefficient function has the same leading small A" singu- 
larities as Cabf(AO) but it vanishes as A" — > oo. It is given by 



C 



(n) 

ABF-sub 



2C^ F (N+l) + C<g> F (N + 2) 



(2.54) 



It is apparent that CABF-sub ("^0 an( ^ ^ABF have the same leading A" 
only introduces subleading N = and A" = —1 singularities. However, limjv 
course, many forms of the subtraction are possible: the particular one given in Eq. (2.54) has been 



*(») 



1 singularities: the subtraction 

(n) 

ABF-sub 



>00 A 



AN) = 0. Of 



chosen as a compromise between the contrasting goals of not changing the small A" behavior and of 



damping strongly enough at large AT. In z space, the subtraction Eq. (2.54) corresponds to damping 
the z — > 1 behavior of the coefficient function through a multiplicative factor (1 — z) 2 . 

In view of combining the small and large N approximations to the coefficient function, one may 
ask what is the expected transition point between the two approximations. In order to answer the 
question, a relevant observation is to note that momentum conservation implies that 7 + (2) = to 
any order in perturbation theory. This in particular implies that all C^ F (N) vanish at N = 2. This 
suggests that N = 2, which is a fixed point for the anomalous dimension, marks the transition between 
the small A^ approximation (not accurate when N >2) and the large N approximation (not accurate 



when A^ < 2). In particular, because the coefficient function Eq. (2.53) is a polynomial in j + (N), it 
vanishes at N = 2 if the anomalous dimension does. 



However, the small N approximation Eq. (2.50) to the anomalous dimension does not respect 
momentum conservation, because it only includes the contribution to 7 + from the leading and next- 
to- leading singularities in A^ = 1 Eq. (2.51 ), and not the full fixed order expression of 7 + Eq. (2.51). 



Momentum conservation can be enforced |47| by adding to C^ F (N) a function / mom (JV), This 
function must not introduce spurious singularities at A" = 1 and it should also be subdominant with 
respect to the large N contributions that we control in C so f t (N). A natural choice appears to be 
/mom(A^) = c/N, with c fixed so that, after subtraction Eq. (2.54), our small A^ coefficient function 
vanishes in A^ = 2. With this choice, the small N approximation of the coefficient function becomes 



C<&(N) = C 



(n) 

ABF-sub 



(AT) 



4' k 



N(N + l)(N + 2)' 



(2.55) 
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Note however that the exact coefficient function does not in general vanish at A = 2, only the 
contribution to it driven by hard radiation from external legs and expressed in terms of the anomalous 
dimension does. Thus, for instance, contributions from subdominant poles in A = 0, —1, —2, . . . will 
in general lead to a non- vanishing contribution to the coefficient function in A = 2. Because we do 
not control such a contribution, we estimate it by allowing the coefficient function to deviate from zero 
at iV = 2, by modifying the value of the constant in the subtraction term of Eq. (2.55). We take this 
deviation from zero to reach as its maximum value 5% of the size of the soft contribution Eq. (2.40) 
! Csoft(2), with either sign; namely, we choose in Eq. (2.55) 



at A 



C AB F- S ub(2)±0.05xC7 soft (2). 



This means that the small A contribution, rather than being completely switched off at A 
small at that point, and gets switched off somewhere in its vicinity. 



(2.56) 



2, is 



Our final result Eq. (2.55) for the small A contribution Cu" (A) to the coefficient function, as well 



as several small A approximations are compared to each other and to the known full result at NLO 
and NNLO in Fig. [2j The full result is shown both in the pointlike approximation (labelled as large 
mt), and for finite mt (labelled as exact): the different small A behavior of the pointlike result, due 
to spurious double poles, is apparent. The small N approximations Eq. ( 2.53| ) (labelled ABF), and its 
subtracted version Eq. (2.54) (labelled ABF-sub) are seen to provide an equally good approximation 
to the exact result in the very small N region where the latter is dominated by its small A poles, but 



only the subtracted version vanishes at large A. The final result Eq. (2.55) after enforcing momentum 



conservation of the anomalous dimension is finally shown (labelled high energy), with an uncertainty 
band obtained by varying the size ofC<2 (2) about zero as discussed above: it coincides with the small 
iV approximation for 1 < N < 1.25, but it is gradually switched off for larger iV until vanishing in the 
vicinity of A ~ 2. 



3 Approximate cross sections up to N 3 LO 



3.1 Parton level results 



We can now construct an approximation to the full coefficient function. Having constructed a large 
A approximation C so ft(A, a s ) Eq. (2.40) and a small A approximation Ch. e . (A, a s ) Eq. (2.55) to the 
coefficient function, in such a way that the small A term does not spoil the large A singularities and 



conversely, we can combine them using Eq. (2.8), which we then expand out in powers of a s according 
to Eq. so that at N fe LO we have 



C (fc) (A) 



(fc) 



(3.1) 



Before turning to the N 3 LO, which is our main result, we first compare the NLO and NNLO 
results found using our procedure to the corresponding exact results. We will use m# = 125 GeV 
and mt = 172.5 GeV t hrou ghout. The comparison is shown in Fig. [3j where our best approximate 



result Capprox(A0 Eq. (3.1) (labelled as approx) is shown along with the large A cj£L(N) (labelled 



as soft) and small A Cj^(A) (labelled as high-energy) terms which contribute to it. As discussed in 



Sect. 2.1 and Sect. 2.2 respectively, the uncertainty on C so ft(A, a s ) is obtained as the spread between 



the two different forms Eq. (2.36) of the large A approximation (green band), while the uncertainty 



on Ch.e. (A, &s) is obtained by varying the size of c£2(2) about zero (blue band). The uncertainty on 



cipp r ox(A) Eq. (3.1) is then obtained as the envelope of these uncertainty bands (red band). In each 



plot we also show the ratio of the approximate result to the exact one. 

It is apparent that the approximate results reproduce the exact one within the uncertainty in the 
full region of real A > 1 at NLO, while at NNLO there is a small disagreement (of about 5%) very close 
to A = 1. Note, however, that in this region what we call "exact" result is not necessarily reliable: 



13 



Higgs partonic gg. Order a s Mellin transform 



Higgs partonic gg. Order a s Mellin transform 




Figure 3. Comparison of the NLO (left) and NNLO (right) exact coefficient functions to various approx- 
imations to it. The large N approximation, corresponding to the band between the softi and soft2 curves 
in Fig. [I] (soft); the small N approximation, corresponding to the high energy curve in Fig. [2] (high energy); 
and the combined small and large N approximation Eq. (3.1 1 (approx). The bottom plot shows the ratio of 
the approximate results to the exact result. Note that at NNLO the "exact" result is in fact the approximate 



construction of Ref. 11 12. 



Higgs partonic gg. Order a s Mellin transform 



2000 




soft 

high energy 
approx 



4 5 
N 



Figure 4. Same as Fig. 3l but at N 3 LO order 



indeed, in the absence of a full NNLO result we are taking as exact the matching of Ref. [IT 12 
of a double expansion in powers of 1 — z and mjj/mt with the exact leading (double) N = 1 pole 
computed in Ref. (To). In particular, the contribution from subleading poles (single pole at N = 1 and 



multiple poles for non-positive integer N) in the "exact" result are not correctly reproduced, while in 
our approximate expression they are partly estimated by varying the size of Cu" (2). 



We now consider our new result for the N 3 LO coefficient function: C, 

(3) 



(3) 



Eq. ( |2.36a > or Eq. ( 2.36b[ ), and (N) is given in Eqs. 



The coefficients in the large N contribution Eqs. (2.36 

coefficient (70,3 > which is unknown: unless stated otherwise, the results are presented with go, 3 



so ft (AO is given either by 
2T54J), g55b in terms of C^ F (N) Eq. \2.53\. 



are collected in Appendix A.l| except the 

0. 



This is a coefficient in the expansion of the constant function g~Q(a s ), related by Eq. (2.37) to the 
function go(oi s ) which appears in the resummed expression Eq. (2.9). The general relation between 
the coefficients go,n and g~o,m is discussed in Appendix A. 1 see in particular Eq. (A. 17): it turns out 
(see Tab. [2] and Ref. |28| ) that the known coefficients go, n are rather larger than go,n- This is also 
the case for Drell-Yan production [28]. For this reason, at third order, where 50,3 = <?o,3 + r 3 ; with 
r% = 114.7, we will take go,3 = (rather than 50,3 = 0) as preferred choice (see also a corresponding 
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Figure 5. Behaviour of the perturbative expansion of C(N, a s ), with a s {m 2 H ) = 0.1126. 



discussion in Ref. 



Coming now to the small N expression Eq. (2.53), the coefficients Cij are 



A.2 



while explicit expressions for are given in Eq. (2.51) 



CaDDrox(iV) is plotted in Fig. ij together with the soft approximation (bounded by the 



The function 
two curves 



collected in Appendix 

,(3) 

y approx 

("3) (q\ 

(N) and (N)) and the high energy approximation (given by C } 

We finally turn to the behaviour of the perturbative expansion of the coefficient function. In Fig. [5] 
we compare the NLO, NNLO and N 3 LO truncations of 



C N 3 LO (iV, a s ) = l + a s CU(N) + a 2 C^(N) + a 3 s Cg » (JNT) 



(3.2) 



We note that at moderately large N > 4 (where we expect our approximation to be very accurate) the 
O(af) contribution is significant, so the convergence of the series is quite slow. On the other hand, 
the saddle point argument of Ref. [29] implies that the dominant contribution at LHC energies comes 
from the region ~ 2, where convergence is much faster, though the N 3 LO contribution is still quite 
large. 



3.2 Hadron level results 

We now discuss the corresponding hadron level quantities. To this purpose, we define the gluon 
channel i^-factors 

K 99 (t, m%) = ll^J^ = 1 + *sK$ + o2*g> + o?.K$ + 0(a% (3.3) 

where a s = a s {m 2 H ), o g g(T,m 2 H ) is the contribution from the gluon channel to the cross section 
Eq. (2.1), which implies that 

<7 (o) (T,m|0 =Ta (m 2 H ,a s )^f gg (T,m 2 H ), (3.4) 

and 

and we use everywhere the NNLO expression of a s , and NNLO parton distributions. 

We then compute the -fT-factors using various approximations for the coefficient function, and 
compare them to each other and, at NLO and NNLO, to the exact result. Specifically, besides our 
preferred approximation Eq. (3.1) shown in Figs. ^H4l we also show results obtained using the soft 
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Figure 6. The NLO, NNLO and N 3 LO contributions to the if -factor Eq. (3.5) in the gluon channel only as 
a function of the collider energy ^fs 1 computed using the various approximations to the coefficient function 
shown in Figs. [3]|4j We also show approximation based on using the TV-soft coefficient function of Fig. [T] 
which at N 3 LO is close to the result of Ref. [28] . 

contribution C^ t to the coefficient function (also shown in Figs. 3-4]), as well as the iV-soft approx- 
imation Eq. (2.35) shown in Fig. [I] and also determined with go,3 = 0. This iV-soft approximation is 
essentially the same as the N 3 LO approximation previously published in Ref. |28|, though here, unlike 
in Ref. 28 , we include the full mt dependence. We use the NNLO NNPDF2.1 [53] set of parton 
distribution function, with a s (m 2 z ) = 0.119. For the time being, we set the scales to their central 
value fi F = /u R = ttih- 

Results are shown in Fig. [gJ where the various contributions to the functions Kg^ are plotted as a 
function of the collider energy yfs. As the energy increases, r becomes smaller and one would expect 
small z effects to become more relevant. Indeed, we observe that the soft approximations deviates from 
the exact results, while the full approximation reproduces well the shape of the cross section for all yfs. 
On the other hand, at low energies the red and green curves (and corresponding uncertainty bands) 
tend to coincide, meaning that the small z contribution has become negligible. In that region, we also 



observe that the bottom edge of the uncertainty band (which is obtained using C so ft 2 Eq. (2.36b)) 



better approximates the exact result than the top (obtained using C ao ft 1 Eq. (2.36a)). Finally, we note 



the iV-soft curve always undershoots the exact result, the more so at higher perturbative orders. This 
agrees with the behaviour of the iV-soft curve in Fig. [TJ 

In Fig. [6] we also show our prediction for the N 3 LO if-factor. The third order term is quite large 
at small energy ^/s ~ 1 TeV, where the soft contribution is dominant, but it remains sizable even for 
y/s ~ 10 TeV, where it is almost half of the NNLO. This slow convergence of the perturbative series 
may make all order resummation mandatory. In the low energy region the uncertainty band on our 
soft approximation is quite narrow, and the TV-soft curve is well below and outside it. At very high 
energies the convergence of the perturbative expansion seems to improve, but very slowly, and the 
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r 



(n) 



exact 



a 



(«) 



approx 



(n) 

TV-soft 



a,K, 



(i) 



a 3 K, 



an 

(2) 

.to 

(3) 
39 



1.328 1.330 ±0.099 1.241 
0.903 0.968 ±0.088 0.815 
0.527 ±0.043 0.353 



Table 1. NLO, NNLO and N 3 LO contributions to the gluon fusion if -factors Eq. (3.3). Our approximation 



is compared to the exact result (when available) and to the TV-soft approximation (close to the results of 
Ref . 28 ) . The uncertainty shown corresponds to the band in Fig. [61 



uncertainty on our prediction increases. 



4 N 3 LO Higgs production at the LHC 



s = 8 TeV, with m H = 125 GeV. In Table [T] 
Eq. computed again with the NNLO NNPDF2.1 [53] PDF set 

0.119, and with our central scale choice fi F = /% = run- Results at NLO, NNLO 



We now concentrate on Higgs production at the LHC at 
we present the iT-factor K gg 
and a s {m 2 z ) 



and N 3 LO are compared to the exact results (when available) as well as the TV-soft approximation 
which, as mentioned in Sect. 3.2, is essentially the same as the approximation published in Ref. [28] . 
The uncertainty on our prediction has been determined as discussed in Sect. [3j Our approximate 
result agrees with the exact result at NLO and NNLO within its stated uncertainty. The iV-soft 
approximation leads to a systematically smaller result, the more so at higher perturbative orders. 

We now turn to our result for the full N 3 LO Higgs production cross section at central scales 



a 



N 3 LO 
approx 



(r,m 2 H ) = a^(r,m 2 H ) 



approx 



(22.58 ± 0.27 + 0.92 • 10~ 2 g ,3) pb, 



(4.1) 



where the error shown is our estimate of the uncertainty in our approximation procedure, and we 
have separated off the contribution from the unknown coefficient g~0 3> discussed in Sect. 3.1 above. 



The LO, NLO and NNLO contributions to the cross section have been computed using the code 
ihixs [54] (corrected for finite top mass effects |6j[T2]) and including all partonic subprocesses, while 
the prediction at N 3 LO only contains the approximate coefficient function for the gluon channel. The 
N 3 LO amounts to a 17% correction to the NNLO prediction cjnnlo = 19.28 pb. This correction 



is somewhat larger than that found in Refs. 55,56 using NNLL resummation, which increases the 
NNLO result by about 8%. Part of the difference is due to the fact that in the NNLL expression 
go(a s ), Eq. (2.9) is included up to order a 2 . Hence, the resummed prediction should be closer to our 
result with go, 3 = 0, which indeed would only lead to a 12% increase over the NNLO. The rest of 
the difference is likely due to the same reasons why our result differs from the iV-soft approximation, 



although the result of Refs. [55 56 differs from the iV-soft construction also by subleading terms of 
collinear origin. 

We now study the dependence of the cross section on variations of the renormalization and fac- 
torization scales, /% and fi F respectively. We first show the scale dependence of the known NLO 
and NNLO cross sections in Fig. [7j We consider both a simultaneous scale variation in all partonic 
subprocesses (black curves), as well as the scale variation for the gluon-gluon subprocess only. The 
renormalization scale dependence of the full result is not much different from that of the gluon contri- 
bution. The factorization scale dependence instead is much stronger for the gluon channel alone than 
for the full result. This cancellation of the factorization scale dependence between partonic subchan- 
nels is a direct consequence of the known structure of the Altarelli-Parisi equations. The factorization 
scale of the full result turns out to be much weaker than its renormalization scale dependence. 
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Higgs hadron-level cross section 



Higgs hadron-level cross section 




Figure 7. Dependence of the NLO and NNLO cross sections on the renormalization scale /i R and fac- 
torization scale ji F . The curves labeled gg (N)NLO are obtained including all channels at (N)LO and the 
gluon-gluon contribution at (N)NLO. 
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Figure 8. Dependence of the N 3 LO cross section on the renormalization scale fi R . 



The scale dependence of our N 3 LO result is displayed in Fig. [8j We only show the renormalization 
scale dependence: the factorization scale dependence of the N 3 LO result will be weaker than that of 
the NNLO, which is already negligible. Also, our N 3 LO result only includes the (dominant) gluon 
contribution, so its factorization scale dependence would be misleadingly large, and canceled by a 
contribution from the quark channels. 

The N 3 LO contribution reduces the renormalization scale dependence of the NNLO QCD result 
from ±10% to ±7% if the scale is varied in the range 0.5 < fXn/rnjj < 2. We also show the prediction 

(3) 

obtained using the soft approximation C N _ soit , with go,3 = 0, i.e. essentially the approximation of 
Ref. [28 1 . The scale dependence of the iV-soft result is similar to that of our own, but its absolute size 
is rather smaller. Note that in Ref. [28] a smaller value of a s (mz) is adopted, which would lead to a 
yet smaller result. 
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5 Conclusions and Outlook 



We have determined an approximate expression for the N 3 LO Higgs production cross section in gluon 
fusion, with finite top mass. We have considered the dominant gluon channel only. We have found 
that at y/s = 8 TeV this correction leads to a 17% increase of the cross section for ran = 125 GeV 
and it noticeably reduces the scale dependence of the NNLO result. 

Our approximation is based on combining information on the large iV and small N singularities 
of the coefficient function, which are determined from resummation, while making sure that they do 
not interfere with each other, so that large iV terms do not introduce spurious small N singularities, 
and conversely. Because large N resummation is known to higher logarithmic accuracy than small N 
resummation (N 3 LL the former, and only leading log the latter) our approximation looses accuracy 
at small N, but fortunately Higgs production in gluon fusion is dominated by large N terms down to 



fairly moderate N values 29 



The results presented here can be used to improve the prediction for standard model Higgs pro- 
duction in gluon fusion, and to the very least they provide an estimate of the impact of higher order 
corrections on the currently known NNLO result which is rather more reliable that the commonly 
used scale variation. A public code is available at 



http://www.ge . inf n. it/~bonvini/higgs/ 



While we have concentrated on the dominant gluon channel the inclusion of other partonic channels 
along the same lines is possible. More interestingly, our approach could also be extended to the 



construction of approximate expressions for rapidity distributions 21,41]. Both are left for future 
work, as well as the construction of a fully resummed result, in which the large and small N terms are 
included to all orders in a s , and the extension to other processes, specifically Drell-Yan production. 
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A Explicit results for the coefficients 



A.l Large iV contributions 

We present some results on Mellin transformation of plus distributions which appear in perturbative 
calculations, defined by 



1 dz[f(z)] + g(z)= ! dzf(z)[g(z)-g(l)}, 
J 

where g{z) is any test function, regular in < z < 1. The distributions 

T,,, fln k (l-z)\ 
V k (z) 



1 - z 
ln fc ln± N 



In 



i 



V k (z) = V k {z) + 



hi 



it i 



75 ln fc (l - z) 



1 - z 



1 - z 



(A.l) 

(A.2a) 
(A.2b) 
(A.2c) 
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can be obtained, respectively, as the fc-th ^-derivative of the generating distributions 

d k 



V,(z) 



V h (z) 



£-1 



5=0 



d k 



The Mellin transforms 



M[f]= I dzz N - l f{z) 



of the generating distributions Eq. (A. 3) are easily computed: 

r(jv) 



M 

One finds [57] 

V k {N) = M[V k {z) 



M 
M 

,-*/2 



(1 

In*" 1 - 

z 



1 _ v (l: + 1 



r(0 
r(0 



i 



r(Jv + o r(i + o 
- 1 



T(N-£/2) 1 

r(iv + e/2) "r(i + 



^ og (iV) = A^ 



1 



where we have defined 



k + 
1 



i=o 
fc 

i=o 



A(0 



fc + 1 

fc+i 

3 

1 



:i) ln fe+1 ^ — 
A 



V o) ( 

r (?)(i) T( fe+1 - J ')(iV,0) - A( fe+1 ^')(l) 



r(0 



(A.3a) 
(A.3b) 
(A.3c) 

(A.4) 

(A.5a) 
(A.5b) 
(A.5c) 



(1) \r(N) A( fe+1 ^') (A) - A^ +1 ~ j) (l)j (A.6a) 



(A.6b) 
(A.6c) 

(A.7a) 
(A.7b) 



T (A, f ) = r(JV - e/2) A (AT + e/2) 

and the superscripts in round brackets in T(A, £) denote derivatives with respect to £. We note that 
all the expressions in Eqs. (A.6) are of comparable complexity. They all share the same behavior at 
large N term by term in the sums, which implies that asymptotically they only differ by constant 
terms. In particular, 

lim \v k {N)-V k {N)} =0, (A.8) 

iV-^oc L J 

while the constants d k defined in Eq. ( 2.38| ) are given by 



d k = lim 

TV-5-oo L 



V k (N)-V^(N) 



k + l 



_ r (fc+i) m 



(A.9) 



The right-hand sides of Eqs. (A.6) and (A.9) are easily computed with the help of the recursion 
relations 



(A. 10a) 
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3=0 
k 



3=0 



k\ 1 

2 



2-i 



T + 



-2)3 



T ( fc -i)(iV, O)^-(iV) 



(A.lOb) 
(A.lOc) 



In particular, Eq. (A. 10c) shows that ipj(N) with j odd never appear in T( fc )(iV,0). More details can 



be found in Ref. 57 



We now consider the coefficient function in the soft limit, which is completely fixed by the coeffi- 



cients b nt k and by the function go(a s ) appearing in Eq. (2.35), which we reproduce here: 

oo n 

C ICS (N,a s ) = g (a s ) exp ]T a?£ 6 n , fc V^(N), 



(A.ll) 



n=l k=0 



The coefficients b n ^ depend only on soft gluon radiation, and therefore do not depend on mn or rrif. 
They can be computed within the effective theory in which the top is integrated out and the top loop 
shrinks to a point (pointlike approximation). On the other hand, go(a s ) depends on run/rut, and 
therefore it also depends on whether the pointlike approximation is used or not. The dependence of 
go(ct s ) on the ratio mn /rnt obviously affects logarithmic terms by interference in C res (iV, a s ) but such 
dependence is under control. 

We now list the explicit coefficients b n ^ for n = 1,2, 3. We omit the scale dependence, which can 
be restored by imposing scale invariance of the hadronic cross section. The order a s coefficients are 



4C A 



71 



'1,0 



0. 



At order a 2 s we have 



b 
62,1 

&2,0 



2,2 — —9 



7T- 



n ^2 2 ^ 
~J C A + 3^/ 



7T" 



101 11 

+ ~3 



7 

-( 
2 



+ ^C2 + ^3 )C 2 A +[ 



(A.12) 

(A.13a) 
(A.13b) 
(A.13c) 



Finally, at order we have 28 
121 



h,3 



1 

1 

7T 3 



7T' 



27 



44 ^2 4 ^ 2 
27^"/ + Y7° An f 



445 11 \ 3 /289 

■"27 +yC2)^+i — 



/ 15503 188 



+ 



11 



648 



() C2-llC3 + yCl)Ci + 
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11 



(A. 14a) 
(A.14b) 

(A.14c) 



C2 



108 
155 



G0 Cl - yC2C3 - 6C5 )C% 



23 

36 C3+ 30 C| ) C ^ 



1 



/mi 1 19 . ., . .. 
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729 27 



(A.14d) 
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n go, n r n 5o,n 

1 4.9374 3.7779 8.7153 

2 10.92 29.18 40.10 

3 unknown 114.7 unknown 

Table 2. Numerical values of go and go at various perturbative orders. 



We now turn to the function 



g (a s ) = 1 + ^2 



Ot s 90,n- 



(A.15) 



n=l 



The first two terms of the expansions are known, and for ran = 125 GeV, mt = 172.5 GeV and nj = 5, 
are given b}j^] 

5 o,i = 8.7153; g , 2 = 40.10, (A.16) 



but go,3 is still unknown. The function go(a s ) is related to go(a s ) by Eq. (2.37), which we rewrite here 

r (*+i)(i) 



using the explicit values of the dk, Eq. (A.9): 



5o (a s ) = g (a s )exp 



n=l k=0 



k + 1 



We find 



where 



g (a s ) = 1 + £a"5o,n, 



n=l 



90, n — g0,n 



(A.17) 



(A.18) 



(A.19) 



and the r n can be read off Eq. (A.17) order by order in o, S ' It is interesting to observe that each r n 
depends on goj with j < n; in particular, does not depend on the unknown coefficient 50,3- The 
numerical values of r n for n = 1,2,3 are given in Table [2] We note that r% is of order 10 2 , which is 
the order of magnitude of a naive estimate of 50,3 on the basis of the known values of <7o,ii 9o,2- 



A. 2 Small contributions 



The coefficients t i 2 of the small N singularity, Eq. (| 2.45[ ) , were expressed in terms of single and double 
integrals over the off-shell gluon virtualities in Refs |10[|43| . For rnn = 125 GeV, mt = 172.5 GeV and 
rif = 5, their numerical values for fi F = m# are given in MS by 



ci,o 

C2,0 
C3,0 



2.28 
4.12 
8.64 



Cl,l 
C2,l 



5.66 
10.54 



(A.20) 



Factorization scale dependence can be easily restored by the substitutions 



ci,o ci n + 4 



2 

In— £ 



C2,0 
Cl,l 



C2,0 + Cl n«F + 



f 2 



ci 1 + 2ci,o4 + 



2 The NLO coefficient po,i has been computed n ume rically using the implementation of Ref. 
The NNLO coefficient go, 2 was computed in Ref. 



12 



of the exact result |6] . 

as an expansion in powers of (urn /mi) 2 ; . We have checked that 



truncating the expansion to order 4 the result is accurate at the per mille level. 
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t 2 e 3 

C3,0 -> C 3 ,0 + C 2 ,o4 + Ci >0 y + 

^2 ^3 

C2.1 -> C2,l + (C 2 ,0 + Ci,i)4 + 3ci, y + y . (A.21) 
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